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We show that for any Hilbert-space dimension, the optimal universal quantum doner can be 
constructed from essentially the same quantum circuit, i.e., we find a universal design for universal 
doners. In the case of infinite dimensions (which includes continuous variable quantum systems) the 
universal doner reduces to an essentially classical device. More generally, we construct a universal 
quantum circuit for distributing qudits in any dimension which acts covariantly under generalized 
displacements and momentum kicks. The behavior of this covariant distributor is controlled by its 
initial state. We show that suitable choices for this initial state yield both universal doners and 
optimized doners for limited alphabets of states whose states are related by generalized phase-space 
displacements. 
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I. INTRODUCTION 



One of the main tasks in quantum information process- 
ing and quantum computing is the distribution of quan- 
tum information encoded in the states of quantum sys- 
tems. Assume a quantum system labelled as 1 is prepared 
in an unknown pure state described by a state vector |^)i 
in an A^-dimensional Hilbert space. The task is to trans- 
fer ■partially the information encoded in system 1 into a 
second system in a covariant way. That is, the fidelity of 
the operation should not depend on the particular choice 
of the input state |^)i. In addition, we want to con- 
trol the amount of information transferred from system 
1 to system 2. One of the simplest examples of such a 
transformation is state-swapping, when the state of sys- 
tem 1 is swapped with the (known) state of system 2. 
In this case the complete information is transferred. An- 
other option is to leave the system 1 in the original state. 
These two operations can be performed with unit fidelity 
irrespective of the input state of system 1. We can also 
consider a case intermediate between these two limiting 
cases, i.e., between no transfer and the complete trans- 
fer of information. One interesting version of this inter- 
mediate transformation involves the copying (cloning) of 
quantum information from system 1 to system 2, where, 
after the transformation, each of the systems 1 and 2 has 
the same reduced state, which is itself as close as possi- 
ble to the original state |\l/)i(4'|. In this case we often 
require that the fidelity of the information transfer does 
not depend on the initial state. 

It is now well known that quantum information can- 
not be exactly copied [|lj. This no-cloning theorem has 
important consequences for the whole field of quantum 
information processing ||]. Nonetheless, considerable 
progress has been made in the development of approx- 
imate quantum doners. Most of the effort has focused 



on two types of these, universal doners that copy all in- 
put states equally well Q , and probabilistic doners that 
copy a known set of states perfectly, but do so with a 
probability which is less than one B . In what follows we 
shall concentrate on universal devices. 

Quantum cloning both illuminates the limits imposed 
by quantum mechanics on the manipulation of quantum 
information and can be useful in applications. It has 
been shown to be useful in improving the performance of 
imperfect quantum detectors |14| and in improving the 
performance of certain quantum computations |Tq]. In 
addition, it has been shown in Ref. [|6[ that quantum 
doners can be used as optimal eavesdropping devices on 
the 6-state cryptographic protocol. We should also note 
that recently an interesting cloning experiment has been 
proposed |17|| . Moreover, two experiments have been in- 
dependently reported [nqjl9| this year in which cloning 
of optical fields has been realized. 

Universal doners can be either symmetric or asymmet- 
ric. In a symmetric doner the quantum information is 
divided equally and the output clones are identical. In 
an asymmetric doner one of the clones receives more of 
the input quantum information than the other. Symmet- 
ric doners were first developed to copy qubits 0,0], but 
have been extended to copy states in spaces of arbitrary 
dimension |§-§[ , and it has been proven that these don- 
ers are optimal P,p|jlC|]. The study of asymmetric don- 
ers also began with the consideration of qubits |ll|-[l^] 
and has been recently extended to systems of arbitrary 
dimension 0]. What we shall do here is to exhibit a 
quantum circuit for symmetric and asymmetric doners 
in arbitrary numbers of dimensions. In order to empha- 
size that what these devices do is distribute quantum 
information, we shall refer to them as quantum infor- 
mation distributors (QID's). The circuit consists of four 
controUed-NOT gates, or rather their generalization to 
N dimensions, and its form is the same for any num- 



ber of dimensions. There are two inputs to this circuit. 
The first is the state which supphes the information to 
be distributed between the two outputs. The second acts 
as a program and determines how the information is dis- 
tributed. The infinite-dimensional version of this circuit 
allows us to describe quantum information distributors 
for continuous variables. 

Let us formulate our problem more exactly. Assume 
the original quantum system is in a pure state 



N-l 



N-1 
n=0 



(1.1) 



At the output of the quantum information distributor 
we would like to have two quantum systems each with a 
state described in a covariant form 



M°"^ = (i-/3^)WiW + Cai 



p^*)^(l-a2)|^) 



N' 



2(^1 + —I2 



(1.2) 



where the real parameters a and (3 quantify the amount 
of information which has been transferred from one sys- 
tem to the other. In particular, if /3 = 0, then no in- 
formation has been transfered from the original system, 
while if a = 0, then all of the information in system 1 
has been transferred to system 2. The parameters a and 
(3 are related (see below). From the covariant form of 
the output density operators it follows that the fidelity 
of the information transfer is input-state independent. 
The terms proportional to 1/N in the density operators 
describe the amount of noise introduced into the systems 
at the output by the information transfer process. 

Our task in this paper is to develop a quantum cir- 
cuit for the universal quantum information distributor 
for arbitrary-dimensional quantum systems. In Sec. O 
we start our discussion with the mathematical formalism 
needed to investigate our problem. Then in Section pn 
we present a quantum network for the universal quantum 
information distributor, while in Section IV we general- 
ize the discussion to continuous variables. Finally, in 
Section ^ we summarize our results. 



11. FROM DISCRETE TO CONTINUOUS 
VARIABLES 



In order to make the discussion self-contained we first 
present a brief review of the formalism describing quan- 
tum states in a finite-dimensional Hilbert space. Here 
we follow the notation introduced in Refs. 12^J21|] (see 
also Ref. [p2|). Let the A^-dimensional Hilbert space be 
spanned by N orthogonal normalized vectors \xk) and 
equivalently by N vectors \pi), k,l = 0, . . . ,N —1, where 
these bases are related by the discrete Fourier transform 






(2.1) 



Without loss of generality, it can be assumed that these 
bases are sets of eigenvectors of non-commuting opera- 
tors X and P: 



X\xk) = k\xk) , P\pi) = l\pi) 



(2.2) 



that is, 



N-l 



X- 


N-l 


k\xk){xk 


p - 


-1. 


i\pi){pi\ 



(2.3) 



/=o 



For instance, we can assume that the operators X and 
P are related to a discrete position and momentum of 
a particle on a ring with a finite number of equidistant 
sites [g3[. Specifically, we can introduce a length scale, 
L, and two operators, the position x and the momentum 
p, such that 



x\xk) = Xk\Xk) , p\pi) = pi\pi) 



where 



Xk = L\ 



Pi 



(2.4) 



(2.5) 



The length, L can, for example, be taken equal to 
y^Ti/cum, where m is the mass and cu the frequency of 
a quantum "harmonic" oscillator within a finite dimen- 
sional Fock space (in what follows we use units such that 
n^ 1). 

The sq uare d absolute values of the scalar product of 
eigenkets ( |2.2| ) do not depend on the indices k, I: 



\{xk\pi)f^l/N, 



(2.6) 



which means that pairs (fc, I) form a discrete phase space 
(i.e., pairs {k,l) represent "points" of the discrete phase 
space) on which a Wigner function can be defined p4[ . 
Next we introduce operators which shift (cyclicly per- 
mute) the basis vectors |23] : 



Rx{n)\xk) — \X(k+n)modN) 
Rp{m)\pi) = \p(l+m)modN) , 



(2.7) 



where the sums of indices are taken modulo N (this sum- 
mation rule is considered throughout this paper, where 
it is clear we will not explicitly write the symbol mod N). 
The operators Rx{n) and Rp{m) can be expressed as 
powers of the operators Rx{l) and -Rp(l), respectively: 



RAn) ^Kil), Rp{m)^R^{l). (2.8) 

In the x-basis these operators can be expressed as 

{xk\Rx{n)\xi) = Sk+n,l 

{xk\Rp{m)\xi) = SkAexp\i—mlj . (2.9) 

Moreover these operators fulfill the Weyl commutation 
relation |2§-|2|] 

Rx{n)Rp{m) — expli — mn]Rp{m)Rx{n) ; (2.f0) 

although they do not commute, they form a representa- 
tion of an Abelian group in a ray space. We can dis- 
place a state in arbitrary order using Rx{n)Rp(rn) or 
Rp{m)Rx{n), the resulting state will be the same — 
the corresponding kets will differ only by an unimpor- 
tant multiplicative factor. We see that the operators 
Rx{n) and Rplm) displace states in the directions x and 
p, respectively. The product Rx{n)Rp{m) acts as a dis- 
placement operator in the discrete phase space (k, I) P9[ . 
These operators can be expressed via the generators of 
translations (shifts) 

Rx{n) = expf-i— nPJ = exp(-ia;„p) 

Rp{m) — expli—mXj ~ eyip{ip„ix) . (2-11) 

We note that the structure of the group associated with 
the operators Rx{n) and Rp{m) is reminiscent of the 
group of phase-space translations (i.e., the Heisenberg 
group) in quantum mechanics |pO|| . 

A general single-particle state in the 2;-basis can be 
expressed as 



N-l 



N-1 



|*)i = ^cfc|xfc)i; ^|cfcp = i. 



(2.12) 



fc=0 



fc=0 



The basis of maximally entangled two-particle states (the 
analogue of the Bell basis for spin-^ particles) can be 
written as llSOll 



|^m«) = -y= ^ exp(^i—mkj\xk)\xi^k-n)raodN) , (2.13) 

where m,n = 0, . . . ,N — 1. We can also rewrite these 
maximally entangled states in the p-basis: 

1 ^""^ 2 

|2™„) = -y= ^ exp(^-i— n?j|p(„,_;)inodJv)b;) •(2-14) 



N 

t — u 

The states \^mn) form an orthonormal basis 
with 



(2.15) 



AT-l 
m,n— 



= 1(g) 1 



(2.16) 



In order to prove the above relations we have used the 
standard relation X]n=o exp[27ri(fc — k')n/N] — N6k,k'- 

It is interesting to note that the whole set of iV^ max- 
imally entangled states |Sm„) can be generated from the 
state |5oo)23 by the action of local unitary operations 
(shifts) of the form 



|^m«)23 = 12 (8'i?2:(n)i?p(m) I ^00)23 , 



(2.17) 



acting just on system 3 in this particular case. 

From the definition of the states |Smn)23 it follows 
that they are simultaneously eigenstates of the operators 
X2-X3 and A + A: 



{X2 — ^3)|^mn)23 = '^|^mn)23 
{P2 + -F3)|^mn)23 = 'm'\^mn)23 



(2.18) 



We easily see that for TV = 2 the above formalism reduces 
to the well-known spin-i particle (qubit) case. 

Now we introduce generalizations of the two-qubit C- 
NOT gate (see also Ref. [Kl| ) . In the case of qubits the 
C-NOT gate is represented by a two-particle operator 
such that if the first (control) particle labelled a is in 
the state |0) nothing "happens" to the state of the sec- 
ond (target) particle labelled b. If, however, the control 
particle is in the state |1) then the state of the target is 
"flipped", i.e., the state |0) is changed into the state |1) 
and vice versa. Formally we can express the action of 
this C-NOT gate as a two-qubit operator of the form 



D,j,= 



\k)a{k\(E)\{m + k)mod2)b{m\ . (2.19) 



k,m—0 

We note that in principle one can introduce an operator 
Dli. defined as 



Dib= J2 \k)a{k\®\{m-k)iaod2)b{m\. (2.20) 

k,m—0 

In the case of qubits these two operators are equal. This 
is not the case when the dimension of the Hilbert space is 
larger than 2 |31|. Let us generalize the above definition 
of the operator D for N > 2. Before doing so, we shall 
simplify our notation. Because we will work mostly in 
the x-basis we shall use the notation \xk) = \k) where it 
may be done so unambiguously. With this in mind we 
now write 



Dr, 



N-1 

E 

k,m—0 



\k)a{k\(g>\{m + k)inodN)b(m\ . (2.21) 



From the definition (2.21) it follows that the operator 
Dab acts on the basis vectors as 



Dab\k)\m) = |fc)|(fc + m)mod7V) , 



(2.22) 



which means that this operator is equal to the condi- 
tional adder [g^,^ and can be performed with the help 
of a simple quantum network as discussed in [p2| . 

If we take into account th e de finition of the shift op- 
erator Rx{n) given by Eq. (2.7) and the definition of 
the position and momentum operators x and p given by 
Eq. (2.8) we can rewrite the operator Dab as: 



N-l 



Dab= Y. \k)a{k\^Ri'Hk)\m)b{m\ (2.23) 

k,m—0 
N-l 



-^XaPb 



fc=0 



and analogously 



N-l 



D\b^ Y. \k)a{k\<»\{m-k)modN)b{m\ (2.24) 

k,m—0 
N-l 



^iXaPb 



fe=0 



where the superscripts a and b indicate on which Hilbert 
space the given operator acts. Now we see that for A^ > 2 
the two operators D and & do differ; they describe con- 
ditional shifts in opposite directions. We see that the 
generalization of the C-NOT operator are the conditional 
shifts. The amount by which the target (in our case par- 
ticle b) is shifted depends on the state of the control par- 
ticle (a). 



A. Continuous limit 

In the TV — > oo limit we have to take special care in 
handling the expressions for the eigenstates of the posi- 
tion and momentum operators p4| . To avoid divergences 
we have to regularize our states by "smearing" them. In 
other words, the eigenstate of the operator x is replaced 
by a squeezed displaced state (see e.g., Ref. [g9[) with 
reduced quadrature fluctuations in the x direction (see 
below). To express these states explicitly we utilize the 
Wigner function representation, which for pure states is 
defined as 



W.^)ix,p) 



1 



/2^ 



dzip^,^{x-)^*(x+)e^P\i2.25) 



where x± = {x ± z/2) and ip (x) is the wave function 
of the state |^), i.e., (p^^.{x) — (a;|^), which in the |a;) 
basis is expressed as 



l*) = 



dxip^^Ax)\x) 



(2.26) 



where we have used the relation {x\y) = \'2tt6{x — y). 
The Wigner function T/F|^)(x,p) is a quasi-probability 
distribution in phase space and is normalized so that 



1 

2^ 



dxdpW {x,p) = 1 



(2.27) 



where dx dp/2'K is the invariant measure in phase space 
(here we have taken ?i = 1). 

With these definitions we can represent a regularized 
version of an eigenstate of the position operator x with 
mean value equal to zero, \xq), as a state described by 
the Gaussian Wigner function: 



^|xo>(2^'P) =2exp(- 



2f 2 

-e ^x 



V) 



(2.28) 



for which the variances of the position and momentum 
operators are (Ax )^ = \e~'^^ and (Ap)^ = \e^^ , respec- 
tively. The state ( ^.28 ) is a minimum uncertainty state, 
i.e., (Ax)(Ap) = i irrespective of the value of the squeez- 
ing parameter ^. For the mean excitation number we 
find the expression n — sinh ^. We see that in the limit 
^ ^ oo, the state described by Wigner function ( 2.28| ) 
is indeed a state with no fluctuations in the x direction 
at the expense of inflnite fluctuations in the p dire ction. 
In other words in the limit ^ — > oo the state ( 2.28 ) is an 
eigenstate of x. 

Analogously a regularized eigenstate of the momentum 
operator is described by the Wigner function 



=2C„2n 



W^^^^{x,p) = 2exp(-e-^^a;^ - e^V) 



(2.29) 



where for the variances of the position and momentum 
operators we find (Ai)^ = ^e^^ and (Ap)^ = \e~'^^, 
respectively. 



Th e wa ve functions corresponding to the states (2.28) 
and (|2.29D read 



2? ^2, 



^„,(x)=2VV/^exp(-l^), 



and 



'P 



|P0> 



(a;) = 2i/4e-«/2exp(- 



e-2«x2 



(2.30) 



(2.31) 



respectively. We denote the corresponding ket vectors 
as |a;o(C)) and |po(C))j where we have explicitly indi- 
cated that these states are regularized versions of two 
specific eigenstates of the position and momentum oper- 
ators. The ^-distribution of the state 09, > ix) is defined 

as usual, i.e., P, , (x) = \w, , (x)P and is normalized to 

unity as J dx P, .{x)/y2/K = 1. This distribution can 



also be obtained from the Wigner function (2.28) via in- 
tegration over p, i.e.. 






1 



'^P^kn>(2^'P) 



which in the large ^ hmit gives Hni^_ 



(2.32) 



'P^.,A^) 



J 



A/27r5(a;), as expected. 

In an analogous way we define a maximally entangled 
two-mode state |Soo(C)) in t^^e continuous hmit. Specifi- 
cally, we define this state in a regularized form for which 
the Wigner function reads pa] 



M^,- 



^oo> 



(xi,pi;a;2,P2) = 4exp<{ -— \(xi - X2f ^ {pi + P2f] 



-2« 



[(xi +a;2)^ + {pi -P2)^] 



r 



This is a Wigner function describing a two-mode 
squeezed vacuum. If we trace over one of the modes, i.e., 
if we perform an in tegra tion over the parameters X2 and 
P2 we obtain from (2.33) a Wigner function of a thermal 
field 

1 f°° 

W^th(2;i,Pi) = ^ / dx2dp2W^-^^^{xi,pi;x2,P2) 



l-|-2n 



exp 



where n = sinh ^ is the mean excitation number in 
the two-mode squeezed vacu um u nder consideration. We 
note that the thermal state (2.34) is a maximally mixed 
state (i.e., with the state with the highest value of the von 
Neumann entropy) for a given mean excitation number. 
This means that the pure state (2.33) is the most en- 
tangled state for a given mean excitation number. From 
this it follows that to create a truly maximally entangled 
state, i.e., the state (2.33) in the limit ^ —^ oo, an infinite 
number of quanta is needed and so infinite e nergy . 

The two-mode wave function of the state (2.33) in the 
^-representation reads 



"Pr. 



'-oo) 



(a::i;a;2) = V2expf- 



=2? 



^2 



-2C 



-xl 



(2.35) 



where x± = {xi ±0:2). In what follows we shall denote 
this regularized version of the maximally entangled state 
in a semi-infinite Hilbert space as |Soo(^))- Now that we 
have laid out the formalism, we can resume our discus- 
sion. 



III. NETWORK FOR QUANTUM 
INFORMATION DISTRIBUTION 



We have shown earlier [p[|q,|36|| that to perform quan- 
tum cloning we need apart from systems 1 and 2, between 



(2.33) 



which the information is shared, an additional quantum 
system 3 which mediates the distribution of the quan- 
tum information. Following this philosophy, we assume 
a quantum information distributor to be a two-particle 
system (2 and 3) each of the same physical type as the 
original system 1. Let us assume that the quantum dis- 
tributor is initially prepared in the most general two- 
particle pure state 



|«')23 = 



JV-1 

E 

m.k—O 



dmk\m)2\k)3 



(3.1) 



In analogy with the quantum computational network 
used in the quantum doner |36[ we assume the QID net- 
work to be 



U123 = D3iDl^Di3Di2 



(3.2) 



with the idea being that the flow of information in the 
quantum distributor, as described by the unitary opera- 
tor (3.2), is governed by the preparation of the distributor 
itself, i.e., by the choice of the st ate ( p^l ). In other words, 
we imagine the transformation ( |3.2| ) a s a universal "pro- 
cessor" or distributor and the state ( p^ ) as "software" 
through which the information flow is controlled. Using 
relation (2.23) we can rewrite the QID transformation as 
(see also Ref. Q) 

U123 = exp[-i(x3 - X2)pi] exp[-i£i(p2 + P3)] ■ (3.3) 

The distribution of information encoded in the original 
particle is performed via a sequence of four conditional 
shifts D. The output state of the three particle system 
after the four controlled shifts are applied is 



\n) 



123 



^31^^1^13^12|*>l|$> 



(3.4) 



The four operators 
|'^)i|TO)2|fc)3 as 



D act on the basis vectors 



D3iDl^Di3Di2\n)i\m)2\k)3 = \{n - m + k)mod N) i |(m + n)modiV)2 \{k + n)modN) 



3 ■ 



(3.5) 



As we shall see, the choice of the state |<5')23 controls the 
flow of the quantum information contained in the state 
|*)i through the QID. 

Before examining this issue, however, it is useful to 
explore the covariance properties of this distributor for 
any choice of |$)23. A device is covariant with respect to 
the transformation U, if application of U to the input, 
i.e., 1^) -^ U\'^) implies that the output density matrix 
representing the pair of outputs, p^°^^>, transforms as |3q] 



p(™t) ^ {/ <^ f7p(°"')c/-i ® U- 



(3.6) 



When examining whether the distributor is covariant 
with respect to transformations of the form Rx{n)Rp{n), 
it is sufficient to confirm this covariant action for 'dis- 
placements' along the x- and p-axis separately, given by 
Rx{n) and Rp{n), respectively. If the state to be dis- 
tributed is Rx{n)\'i')i wc find 

-Rx(n)i|*)i|$)23 -> Rx{n)Mn)2RAn)3\^}i23 , (3.7) 



where |r2)i23 is given by Eq. ( |3.4| ). Similarly, if the input 
state is Rp{n)\'^)i, we have 

i?p(n)i|*)i|$)23 ^ Rp{n)iRp{n)2Rp{~n)3\n)i23 ■ (3.8) 

Combinations of these two 'displacements' act in the nat- 
ural way, so that if we 'translate' the input state by a cer- 
tain amount, the reduced density matrixes of the three 
outputs are translated by the same amount, and if we 
perform a momentum 'translation' on the input state, the 
reduced density matrixes of outputs 1 and 2 are trans- 
lated in momentum by the same amount, while that of 
output 3 has its momentum translated by the opposite 
amount. This implies that this QID is covariant with 
respect to translations and momentum translations, and 
that the fidelities of the output reduced density matrixes 
are unaffected when these transformations are applied to 
the input. 

Having established the covariant action of our distrib- 
utor in arbitrary dimensions for any input state |$)23) 
we now wish to determine how this state affects the flow 
of quantum information in the QID. 
(i) Let us first assume that the QID state |$) is initially 
prepare d in the maximally entangled state |Soo)23 given 
by Eq. (2.14). Taking the original system to be prepared 
in the state (1.1) we find after the QID transformation 



C/l23|*>l|2o 



l*>lls 



00/23 



(3.9) 



that system 1 remains in the original state while the QID 

remains in its initial maximally entangled state |Soo)23- 

This means that even though the three-particle system 

has interacted via four controlled shifts the total state is 

unchanged. 

(ii) Instead, let us assume the QID is initially prepared 

in the product state 



where |po) is an eigenstate of the momentum operator 
with the mean value equal to zero. At the output we 
then find 



C/i23|*)l|2:o)2bo)3 = |*)2|S, 



00/13 



(3.11) 



which means that the information from the system 1 is 
completely transfered to the system 2 while at the output 
the system 1 and 3 are in the maximally entangled state 
|2oo)i3- Note that the output here is a state-swapped 
version of the output in Eq. (^^) . 

Since these two cases realize the two extreme situations 
(no information transfer and complete information trans- 
fer) it is natural to ask what is the action of the QID if it 
is prepared in a linear superposition of the states |Soo)23 
and |xo)2bo)3- Let us take the input state of the QID to 
be 



|$)23 = a|^00)23 + l3\xo)2\Po)3 



(3.12) 



where a and f3 are real parameters. Note, that from the 
normalization condition ('I'l^) = 1 it follows that these 
parameters must fulfill the condition 



z'^^^- 



(3.13) 



When the QID transformatio n is a pplied with the QID 
initially prepared in the state (3.12) the output state be- 



comes 



?7l23|*>l|$)23 = a|*)l|Soo)23 + /?|*)2|Soo)l3 • (3.14) 

Tracing over the systems 2(1) and 3 we find the reduced 
stated for system 1 (2) at the output to be described by 



..(out) 

Pi 



-(out) 
P2 

~(out) 
P3 



2a(3, 



(3^ 



2a(3 



l)p(in) + !_l: 
N ''^ N 



(3.15) 



'-wi""") 



{N - 2a(3) . 
iV2 



|*)23 = |a;o)2bo)3 



(3.10) 



where p^™^ = |^)(^| is the density operator of the origi- 
nal state of system 1, and (p) is the transposed operator. 

Taking into account condition ( 3.13| ) we can directly 
rewrite the last two density operators in the form (pT^). 
This means that QID is the covariant transformation 
which in a controlled way distributes information be- 
tween the two systems. There is a price to pay for this 
covariant information distribution which is reflected by 
the additional noise. 



A. Cloner 



Let us assume that a = f3, i.e., the two outputs ( 3.15 ) 
are equal. In this case QID acts as a u niver sal quantum 
cloner for arbitrary dimensions. From ( 3.12 ) we find the 
initial state of the cloner to be 



use the regularized versions of the states \xo{£,)), |po(C)) 
and |Soo(C)) ^^ introduced in Section^. The input state 
of the QID in the continuous case can then be written as 



mO)23 = a|Soo(0)23 + /3No(e))2bo(e))3 

which in the a;-basis becomes 



1*^23 = ,^,,, , ^, E (l^")2 + \Xra)2) \x^)3 ■ (3.16) 
\/Z(l\ + II „ 



l*(C))23 = ^ / dX2dX3^{x2,X3)\x2)\x3) 



(4.1) 



(4.2) 



With this initial QID state the output of the cloner yields where 
two clones of the form 



.v(out) -(in) 

p) '^sp) 



1 



N 



-1 



J = 1,2. 



The scaling factor s is 



iV + 2 
2{N + 1) 



(3.17) 



(3.18) 



^^{x2,x:i) = ai)\'s.„„){x2,xs.) + Pij\^^^){x2)il}\p„){xz) , (4.3) 

and the Gaussian functions 4'\xn){x 2), V '|pn'> (^3) ; ^'^d 
i l' \ ~.nn') {x2,xz) are defined by Eqs. ( 2.301 ), ( 2.31 ) and 
( ^.35|) , respectively. 

For finite values of squeezing the states |Soo(0)23 ^^nd 
|a^o(0)2|Po(C))3 8'''c not mutually orthogonal, therefore, 
in order to fulfill the normalization condition for the state 



Finally, system 3 of the cloner has a reduced state given |$) the para meter s a and (3 have to fulfill a condition 
by 



analogous to (3.13) 
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a'^+l3^ + 



Aaf3 



v/4 + 2 sinh^ 2^ 



(4.4) 



i.e., this piece of the cloner is left in a state proportional 
to the transposed state of the original quantum system 
plus completely random noise. 



IV. CONTINUOUS LIMIT 

In what follows we make a connection between the dis- 
crete and the continuous case. The role of the controlled 
shifts (NOTs) in the continuous limit is obvious — it is 
a conditional shift down the x-ax is in phase space. Con- 
sequently, the QID operator ( |3.3| ) has a clear meaning in 
the continuous limit. Our goal now is to f ind t he con- 
tinuous analogue of the initial state |$)23 (3.12) of the 



With this initial preparation of the QID, the universal 
(covariant) information distribution of contin uou s vari- 
ables is realized using the network described in (3.2). The 
operator [/123 acts on the basis states |a;i)i|a;2)2|a;3)3 as 



C^i23|a;i)i|a;2)2k3)3 = 1^1)1122)2123)3 , 



(4.5) 



where zi = xi— X2 + x-i, Z2 = Xi -I- 3:2 , and z^ = Xi+x^. 
Assuming that the original system is initially prepared 
in the state 



I*) 
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/27r 



dxi 'ijj{xi)\xi 
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QID. This is rather straightforward: we simply need to the output of the QID becomes 
I 



C/l23|*)l|*(e))23 



(27r)3/2 



dxi dx2 dx3 ij{xi)p{x2, X3) \xi - X2 + 3:3)1 \xi + 2:2)2 \xi + 2:3)3 



(4.7) 



Upon tracing out modes 2 and 3 we obtain from (4/7) the density operator describing the original system at the 
output of the QID: 



Pi 



(out) 



1 



drjdxi dx[ 'iIj{xi)iI;* {x[)IC{xi — x[;r]) \xi + r]){x[ + ry| , 



(27r)3/2 
where the integral kernel /C(xi — x'l; rj) is given by the expression 

K.{xi -x[;r/) ^ — = / dxp-(^— ^i'^ ^'ij P[ 



x-y I . X + V 



(4.8) 



(4.9) 



From the fact that the trace of the density matrix (4.8) is equal to unity we find that the integral kernel (4.£) has to 
fulfill the condition 



dTjJC{0;Tj) = 1 . 

TT J 

The kernel itself can be expressed in the form 

JC(xi\ri) = a^JCi(xi\ri) + /3^/C2(:ri; 77) + a/3/C3(xi; 77) 



(4.10) 



(4.11) 



where we have introduced the notation xi = xi — x'^. Using the explicit expressions for the wave functions describing 
the input state of the QID we find for the kernel functions 



lCi{xi]ri) = e^exp - 
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X < exp 



3e-2« + e2C 
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It is now easy to check that 
1 /"..., . 1 



/2n 
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71^ 



driICi{0;ri) = 
dijICsiO;'!]) = 



d7y/C2(0;7,) = l 



V4 + 2sini?2^ 



(4.15) 



from which it follows that the kernel /C(0; rj) satisfies con- 
dition ( [I.IOD . 

In what follows we utilize the Wigner-function formal- 
ism to analyze the performance of the cloning machine. 
We find a Wigner W{x,p)i''^*' of the output state ( |4.S| ) 
which we express as a convolution of the Wigner function 
Wi {x,p) of the input mode and the Wigner function 
W'^{x,p) of the kernel (f4.1l|): 



w{x,py; 



(out) 
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X W^"^'{x-x',p + p') 



where 



W^ix^p') 
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1^ 



dze'P''IC{z,x'). 



From our definitions it follows that 

W'^{x',p') = a^W'^^{x\p') + l3'^W^^{x' ,p') 



+ a(3W'^^{x,p') 



(4.18) 



From Eq.(4.12) wc find the Wigner function 

W'^'ix',p') 



W'^^ix',p) = e2«exp 



=2? 



— ^{x +p ) 



which in the limit of large squeezing reads 

W'^'{x\p') ^2Tr S{x')S{p'). 



(4.20) 



(4.21) 



The Wigner function W'^^{x' ,p') of the kernel ( 4.13| ) 
reads 



W'^'ix',p') 
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l + 2n 



exp 



(a;'2 -l-p'2) 

2(1-1-277.) 



(4.22) 



where we have used the notation n = sinh" ^ so that 
cosh 2^ =1-1- 2n. We note that this Wigner function 
for large squeezing (i.e. the large 77, limit) is equal to the 

(4 16) Wigner function of a thermal state (2.34) with the mean 
number of excitations equal to 277. ! 

Analogously we can evaluate the explicit expression 
for the Wigner function W'^''{x' ,p') of the kernel (4.15). 

(4.17) This is rather cumbersome, and, since we are interested 
only in the large squeezing limit, we present the corre- 
sponding Wigner function only in this limit 



W'^^ix',p')c^2V2 



exp 



p2C 
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1 iirtc^r I '\ ii, inr- f i- -CiU 1 i which in the large £ limit Can bc formally expresscd as 

where Ir^ (x , p) are the Wigner functions of the kernels ° ^ j t^ 

JC, (j = l,2,3). We can easily check that W^^{x' ,p') ^ 8V2ne-'^ Six') S{p'). (4.24) 

— / dx' dp' W (x',p') = 1 (4-19) Now we can give the explicit expression for the Wigner 

function, Wi {x,p), of the output mode 1 for which 

which is equivalent to the condition (4.10). we find 



<"*^(x,p) 



^M^f "^ {x,p) + ^ I dx' dp' Wi'""^ [x-x\p- p')W^^ (x',p') + 4 V2e-2^a/3M^f "^ {x,p) (4.25) 

27r J 



I 

We note that in the large ^ limit the third term in the Therefore, the output Wigner function is simply a^ times 

the input Wigner function plus 0^ times the convolution 
of the input Wigner function and that of a thermal state. 
Finally we evaluate the fidelity of the QID transforma- 
tion which is defined as 



right-hand side of Eq.(4.25) will vanish due to the factor 
e~^^. Taking into account that in the large squeezing (i.e. 
large n) limit the function W'^^{x,p) i s esse ntially equal 
to a Wigner function of a thermal field ( 2 . 34 ) with a mean 
excitation number of 2n (we will denote this Wigner func- 
t ion a s Wt\x{x,p] 2n)), we can rewrite the Wigner function 
(4.25) at the output of the QID as 



W^^(™*) {x,p)^ a^T^i^'") (x, p) (4.26) 

+ ^ f dx' dp' wl"'\x',p')Wti,{x - x',p~ p'; 2n). 
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In the large squeezing limit the fidelity (4.27) can be ap- 
proximated as 
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r 



The integrals on the right-hand side of this equation 
are less than or equal to 1, so that the entire right-hand 
side goes to zero as 1/tT as tT -^ cx). Therefore, we find 
that in the large squeezing limit, the fidelity of the QID 



is indeed input-state independent, 
^2 = /32. 



and Ti = a while 



A. Universal continuous doner 



If in the limit ^ ^ oo we take a^ = 0^ = ^ then the 
transformations described above describe the symmetric 
doner. In the limit ^ — * oo, the fidelity of the output 
density matrices to the original state is i. This is con- 
sistent with what we expect from the limit N -^ oo limit 
of Eq. ( |3T7| ). 

The fact that in the continuous case the fidelity of the 
copies is ^ is suggestive; it makes one think of a coin 
toss. In fact, one can construct a continuous (univer- 
sal) doner, which is much simpler than the one given 
above but achieves the same fidelity, whose most impor- 
tant component is a fiipping coin H . This doner has two 
inputs, one for the state we wish to clone and one for a 
completely random state (ideally an infinite temperature 
thermal state). What the doner does is to fiip a coin, 
and if the result is heads, the original input state is sent 
to output 1 and the random state to output 2. If instead 
the result is tails, the input state is sent to output 2 and 
the random state to output 1. Assuming that, on aver- 



age, the overlap between the input state and the random 
state is small, this "doner" will clone the input with a 
fidelity of -i. From this we can conclude that for contin- 
uous quantum systems the universal doner is effectively 
a completely classical device. Indeed, one can verify that 
in this limit there is no entanglement between systems 1 
and 2 of the outgoing particles. This is not true in any 
finite dimensional case. Taking this classical distribution 
as a hint, we can see that this type of continuous doner 
is easily generalized to the case of an arbitrary number of 
inputs, Min, and an arbitrary number of outputs, Mout, 
with Mout > Min. In this case the fidelity of cloning is 
just Afin/Mout, which agrees with Werner's result [|lO[ for 
the optimal doner in the infinite dimensional limit. 



V. CONCLUSIONS 

We have shown that for any dimension, the optimal 
universal quantum doner can be constructed from essen- 
tially the same quantum circuit, i.e., what we have is 
a universal design for universal doners. In the case of 
infinite dimensions (which includes continuous variable 
quantum systems) the universal doner reduces to a clas- 
sical device. By contrast, Cerf, et al., have shown that 
if one designs a continuous doner optimized to copy cer- 
tain sets of states, then one can achieve higher fidelities 
than those available to the truly universal doners studied 
here [p^. In particular, they showed that it is possible to 



design a doner that will copy any coherent state with a 
fidelity of 2/3. Their doner also fits within the structure 
of the QID analyzed here if one chooses the initial doner 
state as 



\<P) 



23 



! dx3 exp I 



{-^^)\X2)\X,+X,) 



(5.1) 



It is interesting to note that this doner also produces ap- 
proximate versions of the transpose of coherent states at 
its third output. The transpose of the coherent state \z) 
is |z*). If the input to the doner is \z), then pg " is a 

Gaussian state (that is, (xj/Sg " '\x') is a Gaussian) that 
is concentrated about the point z* in phase space. It is, 
however, more spread out than a coherent state, and its 
fidelity with the actual transposed state is 1/8. 

For continuous systems, these specialized doners will 
be more useful than the universal one. Because of their 
covariance properties, they will clone any two states that 
differ by only a translation in phase space with the same 
fidelity. For example, the fact that the doner in Ref. fS^] 
clones the vacuum with fidelity 2/3 implies that it clones 
all coherent states with the same fidelity. One can eas- 
ily imagine generalizing this result and designing doners 
to optimally clone entire classes of states; if by choos- 
ing the correct |$)23, the doner has been optimized to 
clone a particular state, |^), by covariance it will auto- 
matically be optimal for all states generated from | ^) by 
displacements in phase space. For finite dimensional sys- 
tems, however, universal doners do better than classical 
devices and the simple universal circuit presented here 
shows how they may be constructed. 
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